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O ■ We study the nucleon-nucleon interaction in a chiral constituent quark model 

by using the resonating group method, convenient for treating the interac- 
tion between composite particles. The calculated phase shifts for the S\ and 
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1 Sq channels show the presence of a strong repulsive core due to the com- 



bined effect of the quark interchange and the spin-flavour structure of the 



effective quark-quark interaction. Such a symmetry structure stems from 
the pseudoscalar meson exchange between the quarks and is a consequence 
of the spontaneous breaking of the chiral symmetry. We perform single and 
f— n ■ coupled channel calculations and show the role of coupling of the AA and 

hidden color CC channels on the behaviour of the phase shifts. 
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I. INTRODUCTION 



Many studies have been devoted so far to the understanding of the nucleon-nucleon (NN) 
interaction starting from models which have been considered to be successful in baryon 
spectroscopy. Here we refer to nonrelativisitc quark models in the framework of which 
calculations of scattering phase shifts can be made quantitatively. We can roughly divide 
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these models into three categories. In the first category we consider models based on one- 
gluon exchange (OGE) between quarks. They explain the short-range repulsion of the 
NN potential as due to the chromomagnetic spin-spin part of OGE, combined with quark 
interchanges between 3g clusters (for a review see e.g. []l]-§]])- in addition, the long-range 
part is obtained from the one-pion exchange potential acting directly between two nucleons 
and the medium-range part is introduced phenomenologically as a local central potential || . 

There is a second category, of hybrid models |f-§|, where in addition to OGE, quarks 
belonging to different clusters interact also via pseudoscalar and scalar meson exchange. In 
these hybrid models the short-range repulsion is still attributed to OGE and the middle- 
and long-range attraction is due to meson exchanges between quarks. 

Here we employ a model of the third category where the quark-quark interaction, besides 
the confinement, is due entirely to meson exchanges between quarks. This is the chiral 
constituent quark model proposed in Ref. and parametrized in a nonrelativistic version 
in Refs. flT0|,pTf. There are also semirelativistic versions available, see e. g. ||12|. For the 



present status of the model we refer the reader to Ref. [pi . 

The origin of the model |9|-|r3|] is thought to lie in the spontaneous breaking of chiral 
symmetry in QCD which implies the existence of constituent quarks with a dynamical mass 
and Goldstone bosons (pseudoscalar mesons). According to the two-scale picture of Manohar 



and Georgi [|14| at a distance beyond that of spontaneous chiral symmetry breaking, but 
within that of the confinement scale, the appropriate degrees of freedom should be the 
constituent quarks and the chiral meson fields. If a quark-pseudoscalar meson coupling is 
assumed, in a nonrelativistic limit one obtains a quark- meson vertex proportional to a ■ q X F 
with a the Pauli matrices, q the momentum of the meson and A the Gell-Mann flavor (F) 
matrices. This generates a pseudoscalar meson exchange interaction between quarks which 
is spin and flavor dependent. 

In the following this interaction is referred to as a Goldstone boson exchange (GBE) 
interaction. In the coordinate space the corresponding potential contains two terms. One is 
a repulsive Yukawa potential tail and the other is an attractive contact 5-interaction. When 



regularized ||10| , |11| the latter generates the short range part of the quark-quark interaction. 
The short-range part dominates over the Yukawa part in the description of baryon spectra 
leading to a correct order of positive and negative parity states |L5] . This interaction contains 
the main ingredients required in the calculation of the NN potential, and it is thus natural 
to study the NN problem within this model. In addition, the two-pion exchange interaction 
between constituent quarks reinforces the effect of the flavor-spin part of the one-meson 
exchange and also provides a contribution of a cr-type scalar meson exchange [16[ required 



to describe the middle-range attraction. 

The spin-flavor symmetry structure of the model []9Hl3f , which is essential in describing 
the light baryon spectrum is getting support from the phenomenological analysis of L = 1 
negative parity resonances \ Tf\ . Also 1/N C QCD studies |18j have a consistent interpretation 



in a constituent quark model with pseudoscalar meson exchange interaction. The spon- 
taneous chiral symmetry breaking is responsible for the absence of parity doubling in low 
energy hadron spectrum. In particular it explains the splitting between the negative parity 
state iV*(1535) and the nucleon iV(939). The quark models, as e. g. the OGE model, which 
explicitly breaks chiral symmetry, fail to reproduce the iV*(1535) — iV(939) splitting. Recent 
lattice calculations, which take into account the chiral symmetry of QCD |T9[, were able to 
reproduce the above N* — iV splitting. This brings new substantial support to the model 

a 



This work is a natural extension of the previous studies [20-22]. Rcf. [20| was rather 



exploratory about the role of a spin-flavor dependent interaction in giving rise to a repulsive 



core. Within the parametrization jnj of the GBE model it was found that at zero-separation 
between two 3q clusters the height of the repulsive core is 0.830 GeV and 1.356 GeV in the 
3 S\ and 1 5'o channels respectively. The spin-flavor symmetry and the parametrization ]nj of 
the GBE model favours the | [42] [51] FS ) state which becomes highly dominant. In a better 
basis [flTJ , obtained from single-particle molecular type states, instead of cluster model states, 
the situation is similar, the repulsion being reduced by about 200 MeV in the 3 S"i channel 
and by about 400 MeV in the 1 S'o channel. This is natural because the molecular orbital 



basis gives a lower bound of the expectation value of the Hamiltonian in the six-quark basis. 



In Ref. |22|] an adiabatic nucleon-nucleon potential was calculated based on the model [jTD[ . 
It was found that none of the bases, cluster or molecular, leads to an attractive pocket. An 
attraction was simulated by introducing a a-meson exchange of a similar analytic structure 
between quarks, with that of the pseudoscalar meson exchange. 

Here, instead of [10|, we use the chiral constituent quark model version of Ref. |TT| 



where the GBE interaction is parametrized in a more realistic way. The adiabatic potential 
calculated |22f in the Born-Oppenheimer approximation with this version posesses a small 



attractive pocket, in contrast to that resulting from model [Hj (see Ref. [f2l|j). 

The present study is based on a dynamical approach to the NN interaction, namely 
the resonating group method (RGM) [f23| -|2"5|| , which allows the calculation of both bound 
states and phase shifts. This method has already been used in NN studies with models of 
categories 1) and 2) mentioned above. So far it has been always applied to nonrelativistic 
models, where the wave function of the nucleon can be approximated by an s 3 configuration. 

In the next section we shortly review the Hamiltonian model [JO]]. In Sec. Ill we describe 
the main steps of the resonating group method for bound and scattering states. The 6q basis 
formed of NN, AA and CC (hidden color) states is introduced in subsection III C. In Sec. 
IV we derive the matrix elements required by the RGM method for the typical spin-flavor 
structure of the GBE model. In Sec. V we present the results for the phase shifts in the 
3 S\ and 1 Sq channels and discuss the role of the coupled A A and CC channels on the NN 
phase shifts. The last section is devoted to conclusions. 

II. THE MODEL 



The GBE Hamiltonian considered below has the form [ |TT|j 

2 

tf = £™* + £^r-*G + £ Voonfa) + £ V x ( rij ) , (2.1) 

where Kg is the kinetic energy of the center of mass. The linear confining interaction is 



Vconfi^) = -|a,! • \ c 3 (Cr i3 + V ) , (2.2) 

and the spin-spin component of the GBE interaction in its SUp(3) form is 

Vxfa) = {£ K(r^)Af Xf + j^ V K ( rij )XfXf + ^(r^AfAj + |^(r«)}d5 • oj . (2.3) 



The interaction (|2.3j) contains j = tt,K, 7], and 7/ meson exchange terms and K,(ry) is given 



as the sum of two distinct contributions: a Yukawa-type potential containing the mass of 
the exchanged meson and a short-range contribution of opposite sign, the role of which is 
crucial in baryon spectroscopy. For a given meson 7, the exchange potential is 



a 2 i P -^r P -A 7 r 

V,(r) = |L— !— {^ Af- } , (2.4) 

47T Izrriimj r ' r 



where A 7 = Aq + k/i 7 . For a system of u and d quarks only, as is the case here, the K 



exchange does not contribute. In the calculations below we use the parameters of Ref. []TT 
These are 

m u4 = 340 MeV , C = 0.77 fnT 2 , 

fx n = 139 MeV , fi v = 547 MeV , /v = 958 MeV > 
22 2 

% = % = 1.24, 9 -^ = 2.7652, 

47T 47T 47T 

A = 5.82 fm^ 1 , k = 1.34 , Vb = -112 MeV . (2.5) 

As already mentioned before, the reason of using the parametrization [11], instead of 



T^], as in the previous work p0|-p2|, is that it is more realistic. Its volume integral, i. e. its 



Fourier transform at q = 0, vanishes, consistently with the quark-pseudoscalar meson vertex 
proportional to a ■ q \ F . In addition, this interaction does not enhance the quark-quark 
matrix elements containing lp relative motion, as it is the case with the parametrization 
flO|j . This point has been raised in Ref. ||26| . 

At this stage we wish to stress that the above parametrization gives a good description 
of baryon spectra. We do not change any parameter obtained from the fit JTTJ. Such a 



parametrization is, of course, only effective. However, irrespective of the parametrization, 



the flavor-spin symmetry is essential in this model. There are also semirelativistic versions 



of the GBE model, as for example |T2| but the application of the RGM techniques to 



semirelativistic six-quark Hamiltonians is certainly much more involved. 



III. THE RESONATING GROUP METHOD 



The resonating group method |23| is one of the well established methods used to study 
the interaction between two composite systems. It allows to calculate bound states energies 
and scattering phase shifts. It has been first applied to nuclear physics in the study of the 
nucleus- nucleus interaction [[24|,[25|]. Its application to baryon-baryon systems was initiated 



by Oka and Yazaki |27|| . In a baryon-baryon system, where each baryon is a 3q cluster, 
it takes explicitly into account the quark interchange between the two interacting baryons. 
This comes from the assumption that the total wave function can be written as 

^ = Y.A[§pXp{Rab)\ , (3.1) 

& 

where (3 is a specific channel (here j3 = NN, A A or CC), A is an antisymmetrization 
operator defined below, $^ contains the product of internal wave functions of the interacting 
baryons and Xp{Rab) is the wave function of the relative motion in the channel f3, depending 
on the relative coordinate Rab between clusters A and B. 

The internal wave function of each cluster has orbital, flavor, spin and color parts. In 
$/3 the flavor and spin are combined to give a definite total spin S and isospin / so that one 
has 



MU)MOb)\ si , (3-2) 



$/3 

where £4 = (^1,^2) and £g = (£3, £4) are the internal coordinates of the clusters A and B: 



^ = ^-^2, £ 3 = f 4 - r 5 , 

p _fi + f 2 -2r z -> _r 4 + r 5 - 2f 6 

£2 - 2 ' ^ 4 " 2 ' 

5 n + r 2 + r 3 B r 4 + r 5 + rQ 

tiA = 7y ) -Kb = 5 • [o-o) 
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The functions &(£%), i = A,B are supposed to be known (see later). They are totally 
antisymmetric 3g states in orbital, spin, flavor and color space. The color part is a [l 3 ] 
singlet for N and A states and an octet for C states. Usually the color part of a 3q state is 
not written explicitly. The same statement remains valid for the 6q state which is a [222] c 
singlet in any channel. 

The antisymmetrization operator A is defined by 

^ = 1-EE^, (3-4) 

i=l j=4 

where P^ is the permutation operator of the quarks % and j belonging to clusters .4(1,2,3) 
and 5(4,5,6) respectively. It acts in the orbital, flavor, spin and color space, so it can be 
written as P tj = P^P/jP^P^ where 

P& = \>l • >i + \ , P? 3 = \^-^ + \, P§ = \>$-Xj + \, (3-5) 

with \{ the Gell-Mann matrices of SUp(3) (SUc(3)) and <r$ the Pauli matrices. 

Let us first consider the one channel case. From the variational principle one can obtain 
the equation determining the relative wave function x(Rab) 

J <P + (£a)<P + (£b)(H - E)Am A ) ( p(£ B ) X (RAB)}d 3 Ud^B = , (3.6) 

where H is the Hamiltonian of the six-quark system. As usually (see e. g. Ref. 0) we 
introduce the Hamiltonian kernel 

H(R',R) = j ' <f> + (U)<f> + (i B )5(R! - R A B)HA{<P{U)<P{£B)KR-RAB)}d?UdHBd z R A B 

= H (d) (R)5(R - B!) - H {ex) (R', R) , (3.7) 

and the normalization kernel 

Af(R', R) = j <j)+{U)(j) + {£ B )8{Ri - Rab)A[(/>(U)<I>(Zb)S(R - RAB)]d 3 Ud^ B d 3 R A B 

= A/" (d) (R)5(R - RJ) - N {ex) (R f , R) . (3.8) 

The direct term of the Hamiltonian kernel, 7v (R), consists of the relative kinetic, the 
relative potential and the internal energies : 



nW(R) = -^k + v}$(R)+H int , (3.9) 

where /i = 3m/2 is the reduced mass of the clusters A and B. Then Eq. ( |3.6| ) can be written 
as 

f C(R f ,R) X (R)d 3 R = , (3.10) 

where C(R',R) = H(R',R) - EJ\f{R',R). This is the RGM equation. Using Q one can 
write 

£{R>, R) = [-^ + V}£{R) - E rel ]5(R - R!) - [H iex) {R>, R) - EM {ex \R>, R)\ , (3.11) 

where E Te \ = E — Hi nt is the energy of the relative motion. There are two important steps 
in solving this equation. One is to calculate the Hamiltonian kernel ( |3.7|) by reducing the 
six-body matrix elements to two-body matrix elements. This is discussed in Sec IV. Another 
step is the discretisation of the RGM equation. It is important both for bound and scattering 
states. The discretisation has been performed by using the method of Ref. 

A. Bound states 

Here we briefly describe the discretisation procedure directly applicable to bound states. 
According to Ref. [ES|, the relative wave function x(i?) has been expanded over a finite 



number of Gaussians Xi centered at Ri (i = 1,2,..., A) where Ri are points, here equally 
spaced, between the origin and some value of R depending on the range of the interaction. 
The expansion is 

N 

X (R) = Y, C ^(R) , (3-12) 

with 



Xi(R) = g(R-Ri, y/2/3b) = (^) 3/ V3^-^ . (3.13) 

If g(r, b) is the normalized Gaussian wave function of a quark, given by 



9(r, b) = {\f /A e-^ , (3.14) 

7T0 Z 



from the Jacobi transformations Q3.3|) it follows that the relative wave function is expanded 



in terms of the Gaussians (|3.13|) with the size parameter J2/3b. This method can be applied 



straightforwardly to the bound state problem. The modification necessary for treating the 
scattering problem will be explained later in the next subsection. The binding energy E and 
the expansion coefficients C{ are given by the eigenvalues and eigenvectors of the following 
equation : 

N N 

^HijC^E^NijCj, (3.15) 

3=1 3=1 



where N is the number of Gaussians considered in Q3.12 ). The matrices 



Ha = j <P + {£a)<P + {£b)x(Rab - Ri)H(l - A'MZaMZbMRab - R,)}d 3 Ud^ B d 3 R A B , 

(3.16) 

and 

Nij = J 4> + (£a)<P + (£b)x(Rab - Ri)(l - M)[4>{£aWb)x{Rab - R 3 )]d 3 Ud 3 ZBd 3 R A B 

(3.17) 

are obtained from ( |3.7| ) and ( |3.8| ) respectively. By including the center of mass coordinate 
(Ra + -R_b)/2 and transforming back to r« (i = 1, ...,6) we get the following formulas 

By = / n <f> + (n - f ) n <^ + ^HAifim - f ) n *&, + ^n...^ , 

J k=l Z k'=4 Z 1=1 Z l'=4 Z 

(3.18) 



and 



p 6 p 3 p 6 p 



% = / n <p + (n - f ) n <f> + (n' + %)A[ii ^(ri - tk) n m + ^)]<*v..rfv 6 , 

J fc=l Z fc'=4 Z 1=1 Z i'=4 Z 

(3.19) 



with 0(f) = g(f,b) given by ( |3.14 ). These forms are much easier to handle in actual 



calculations. They allow to reduce the 6q matrix elements to two-body matrix elements. 

9 



Moreover the distances Ri play now the role of a generator coordinate {4| and lead to a 
better understanding of the relation between the resonating group method and the generator 



coordinate method 28 



B. Scattering states 

For scattering states the expansion (|3.12|) holds up to a finite distance R = R c , depending 
on the range of the interaction. Beyond R c , x{R) becomes the usual combination of Hankel 
functions containing the S-matrix. Because practical calculations of both bound state and 
scattering states are done in terms of partial waves, we first give the partial wave expansion 
of Eq. (|3 . 1 2j ) in terms of locally peaked wave functions with a definite angular momentum 



I and projection m: 

N 



Xi m (R) = E Q l) X ^(R)Y lm (R) , (3.20) 

with the explicit form of xi given by 

* ?)(i?) = An ^ )3/4e ~^ R2+Rl) ^W m) ' (3 - 21) 



where %i is the modified spherical Bessel function [Ea|. When we treat the scattering problem, 
the form (|3.21|) holds up to R < R c only. In fact in this case the relative wave function is 
expanded in terms of x ' as 



x {l \R) = T.crr:\R), (3.22) 



i=i 



where 



x!\r) = 4 ) x!\r), (r<r c ) 

X?(R) = h\-\kR) + Sl l) h^\kR) , (R > R c ) (3.23) 



with x\ (R) defined by Eq. ( |3.21| ). Here k is the wave number k = J2fiE re i and h\ 
and h\ are spherical Hankel functions |2P[ . The coefficients a\ and S\ are determined 
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from the continuity of xi an d its derivative at R = R c . The coefficients C\ of Q3.20 ) are 

N 

normalized such that y^C,; = 1. Then the S'-matrix is given in terms of the coefficients 

C/^ as 

N 

s (i) = j2cPsi l) . (3.24) 

The method of determining the expansion coefficients is described in detail by Oka and 



Yazaki |27j 



C. Coupled channels 

Here we consider more than one channel. In this case, based on Eq. (|3.1| ), the RGM 
equation becomes a system of coupled channel equations for xp 

J2 fc af3 (R',R) Xf 3(R)d 3 R = £ f[H a p(R',R) - EM a p(R',R)] X p(R)d 3 R = . (3.25) 
p J p J 

Usually the normalisation kernel M a p is not diagonal because of the antisymmetrisation. 



For a given SI sector one can establish which are the Qq states of ( |3.2| ) allowed by the Pauli 
principle |3(J. Here we consider the I = partial waves i. e. we study the 3 S± and 1 Sq 



phase shifts. In this case, according to JHJ], the 6q allowed states are NN, A A and CC. 
The NN and A A states are easy to define directly from Eq. (|3.1|). For CC states we adopt 
the definition of Ref. [|3l] which is more appropriate for RGM calculation. This CC state of 



six quarks allows some "color polarisation" of the 6q system in the interaction region. It is 
defined in the following way 

\CC) = a\NN) + /3|AA) + 1 A STC \ AA) , (3.26) 

with 

-^=^ 1 -EE^]. ( 3 - 27 ) 

iU i=l j=4 

where P^,P/j and P^ are the exchange operators in the spin, isospin and color space respec- 
tively defined by ( |3.5| ). From the orthonormality conditions (CC\CC) = 1, (CC|iViV) = 
and (CC|AA) = one can determine the coefficients a, (3 and 7 so that 

11 



\CC) = -%NN) + I|AA> - ^Astc\AA) . (3.28) 

o 3 4 

The important feature in the definition of the CC-state is that the eigenvalue of the color 
SU(3) Casimir operator is 12 for each 3q cluster. This tells us that C is a color octet 
state and thus explains why we call the CC-state a hidden color state. Note that at zero 
separation between quarks (shell model basis) the CC state above is the same as that 
introduced by Harvey. The two differ only at finite separation distances. To see the identity 
with Harvey's CC state |30| at zero separation one can combine it with the A^A" and AA 
states as defined by Eq. (|3.1|) to get symmetry states of the form | [/]irs[222]c; (Jfsc) where g 
is the representation resulting from the inner product of [/]fs an d [222]c which is conjugate 
with the symmetry g of an orbital state such as to produce a totally antisymmetric 6g state. 



Comparing Table 3 of Ref. [31] with that of Harvey's [[3(J Table 1 one can see that the 
coefficients of this basis transformation are identical which proves the identity of the hidden 
color state ( |3.28| ) with that of Harvey at R = 0. Note that Harvey's definition J5C] of CC is 



more appropriate for generator coordinate method than for RGM calculations. 

IV. SIX-BODY MATRIX ELEMENTS 

The method to compute the six-body matrix elements is explained in some detail in the 
appendix. In Tables I & II we give the results for diagonal and off-diagonal matrix elements 
of the channels NN, AA and CC to be used in coupled channel calculations of the 3 Si and 
1 S'o phase shifts respectively. Although we apply the SU(3) version of the GBE model the 
matrix elements of o~i ■ Oj r% ■ Tj and a- L ■ Oj Ti ■ rjP 3 g c needed in SU(2) calculations are also 
indicated. In fact they are used in calculating the expectation value of ctj ■ Oj A® ■ A® by 
subtracting them from <7j ■ crj \[ ■ AJ because there is no K meson exchange. Moreover the 
values we found can be considered as a validity test of our method because they are in full 
agreement with Table 1 of Ref. 
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V. NUMERICAL RESULTS 

We perform RGM calculation as described above for NN, NN + AA and NN + AA + CC 
channels. In all cases the size parameter of the Gaussian ( [3.14|) is fixed at b = 0.44 fm by 
the stability condition (see for example Ref. I]) 

^<0|#|0>=O, (5.1) 



where is a variational solution of the Hamiltonian fl2.1|) for a ground state 3q system. This 



solution is fully symmetric in the orbital space and is chosen to be of the form 

<P=f[g(n,b) , (5.2) 

4 = 1 



with g{fi, b) of (|3~T4T). 



Either if we take one, two or three channels namely NN, NN + AA or A^ A" + AA + CC 
we found that a number of 15 Gaussians in the expansion ( |3.12| ) is large enough to obtain 
convergence. In all cases the result is stable at the matching radius R c = 4.5 fm. In Figs. 1 
& 2 we show the phase shifts as a function of the relative momentum k obtained from one, 
two and three coupled channels. One can see that the addition to A^ A" of the AA channel 
alone or of both AA and CC channels brings a very small change in the 3 S\ and 1 S'o phase 
shifts below 2.5 fm , making the repulsion slightly weaker. The CC channel brings slightly 
more repulsion than the AA channel. In fact the role of CC channels is expected to increase 
for larger values of k, or alternatively smaller separation distances between nucleons, where 
they could bring an important contribution. Of course, the contribution of the CC channels 
to the A^ A" phase shifts vanishes at larger separations because of their color structure. The 
conclusion regarding the minor contribution of AA and CC channels to the phase shifts 
below 2.5 fm -1 is similar for results based on the OGE model (see for example |31|). Thus 
for I = waves it is good enough to perform one channel calculations in the lab energy 
interval 0-350 MeV. 

We recall that the pseudoscalar exchange interaction ( |2.4| ) contains both a short range 
part, responsible for the repulsion, and a long range Yukawa- type potential which brings 

13 



attraction in the NN potential. In order to see the difference in the amount of repulsion 
induced by the GBE and that induced by the OGE interaction we repeated the one channel 
(NN) calculations above by removing the Yukawa- type part. We compared the resulting 
phase shifts with those of Fig. 2 of Ref. |H|] obtained with an OGE interaction parametrized 
such as to satisfy the stability condition ( |5.1| ). We found that in the GBE model the repulsion 
is much stronger and corresponds to a hard core radius Tq BE = 0.68 fm (versus Tq GE = 0.30 
fm) in the 3 S± and r = 0.81 fm (versus r = 0.35 fm) in the 1 So partial waves. The 
radius r$ was extracted from the phase shifts at small k, which is approximately given by 
5 = —kr$. One can also see that the repulsion induced by the GBE interaction in the Z S\ 
partial wave is weaker than that induced in the 1 Sq partial wave. This is consistent with 
our previous result |22| where we found that the height of the repulsive core is lower for 3 S\ 
than for 1 Sq, as mentioned in the introduction. Thus the OGE model gives less replusion 
than the GBE model. In Ref. |33[ the stronger replusion induced by the GBE interaction is 



viewed as a welcome feature in correctly describing the phase shifts above Ei a ^ = 350 MeV. 

A note of caution is required regarding the removal of the long-range Yukawa part of 
the interaction (|2.4j) with the parametrization ( |2.5| ) which contains a rather large coupling 
constant g 2 > /(47r) = 2.7652. The r\ -meson exchange is responsible for describing correctly 
the A— N splitting. If the long-range Yukawa part is removed, the model fails to describe this 
splitting because the contribution coming from the second term of ( [2.4|) for j — r\ becomes 
too large in a 3q system in the parametrization ( |2.5|) . We recall that the contribution to iV of 
the short-range rj -meson exchange part is proportional to a factor of 2 and the contribution 
to A to a factor -2 || , which brings A too low and N too high if the Yukawa part is removed. 
In these circumstances two or three coupled channel calculations become meaningless. 

It is also interesting to see the behaviour of the relative wave function x' =0 of Eq. ( |3.22| ) 
at short distances. Instead of \ = ° ft is more appropriate 0] to introduce a renormalized 
wave function as 

X l =°(R) = £ [dti [N l f a °(R,R')] 1/2 X l A R ') - (5-3) 

p J 
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where the quantity to be integrated contains the I = component of the norm N. In 
Fig. 3 we show results for the above function for the 3 Si wave at k — 1 fm _1 both for 
the one and the three channel cases. One can see that for R < 1 f m the two functions are 
entirely different, in the three channel case a node being present. If the renormalization 
was made with the norm N instead of its square, as in Eq. ( |5.3j ), no node would have been 
present. The existence of a node is related to the presence of the [42] o configuration in 
the wave function (see e.g. [20|). Here, whenever it appears, it is due to the cancellation of 



the positive and negative components of the wave function, but the lack of a node does not 
exclude a repulsive potential. In a renormalized wave function the amplitudes of positive 
and negative components change their values depending on the multiplicative factor iV or 
N 1 ' 2 so the node could appear in one renormalization definition but not in the other. On 
the other hand, as discussed above, the phase shift changes insignificantly when one goes 
from one channel to three channels, and this can also be seen in the asymptotic form of 
the wave function beyond R = 1 fm, although in the overlap region the two functions are 
entirely different. The above behaviour of the wave function is very similar to that found in 
Ref. |5tJ where no long-range part is present in the schematic quark-quark potential due to 
pion exchange. 

In Fig. 4 we represent the 3 S\ and 1 S'o phase shifts of Figs. 1 & 2 in the one channel case 
(NN) again with the Yukawa part included, but this time as a function of Ei a f, = 2h 2 k 2 /3m 
with m = m Ut d of (|2.5|) . This is to show that in the GBE model the two phase shifts are 
very near each other, with 5( 3 Si) slightly lower than <5( 1 S'o). Contrary, in OGE calculations 
as example those of Fig. 2 of Ref. |U| one obtaines S( 3 Si) > S^Sq). In calculations based 



on the OGE model the difference between the two phase shifts is reduced by the addition 
of a scalar potential acting at a nucleon level with a larger attractive strength in the So 
channel than in the 3 Si channel ||. 

A major difference between the GBE 5( 3 Si) and 5^ So) is expected to appear after the 
inclusion of a quark-quark tensor force [|34|. This will modify only the 3 S± phase shift. 
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VI. CONCLUSIONS 



This work is a further important step in our previous studies PTLE3 of the A^ A" prob 



lem. We consider the two interacting nucleons as a 6q system described by a Hamiltonian 
confining a linear confinement plus a pseudoscalar (meson) exchange interaction between 
quarks. 

Previously we derived an A^A" potential in an adiabatic approximation. The present 
study is based on a dynamical approach of the A^A" interaction, namely the resonating 
group method. We perform one, two and three coupled channel calculations for the 3 S*i and 
1 5'o phase shifts for laboratory energies up to about 350 MeV. 

Our conclusions are : 

1. The phase shifts present a behaviour typical for strongly replusive potentials. We find 
that this repulsion, which is induced by pseudoscalar meson exchange is stronger than 
that produced by the OGE interaction. 

2. In the 1 S'o partial wave the repulsion is stronger than in S% partial wave as our previous 
studies suggested. 

3. Our results prove that in the laboratory energy interval 0-350 MeV the one channel 
approximation is entirely satisfactory. 

Finally in future calculation, in order to describe the 3 S± phase shift the tensor force is 
compulsory and this is our following major step. 



Acknowledgements. We are most grateful to Kiyotaka Shimizu for help in understanding 
the resonating group method techniques and for constructive criticism in preparing the 
manuscript. 
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The method to compute the six-body matrix elements is explained here using the example 
of S = 1, 1 = case. 

We know that for the nucleon, the spin-flavor wavefunction is given by 

^= * lyy' + xV], (o.i) 

where x an d 4> are the spin and flavor parts respectively. For the spin parts we have 

Xi/2 = ^(TIT - ITT) , 
X-i/2 = 7|(T1! ~ i-TI) , 

Xi/2 = 4?(UT + ITT -2 TU) , 
xi 1/2 = ^(TU + ITI -2 UT) , (0.2) 

and similarly for the flavor parts with | replaced by u and J. replaced by d. Then for 
/3 = NN, the Eq. (O) becomes 



*w = ^Ectt^&K(l)<(l) H-X S A 1 (1)^ 1 (1)]K(2)<(2) H- X A 2 (2)<(2)] , (0.3) 

where 5 and I are the spin and isospin of the NN system. x(z) and <p(i) are the spin and 
flavor parts of the i nucleon. For S = S z = 1 and I = I z = 0, after inserting the values of 
the corresponding Clebsch-Gordan coefficients we have 

®n\ = ^{[Xx/ 2 (l)0i/ 2 (l) + X^ l)^/ 2 (l)M/ 2 (2)C 1/2 (2) + X A /2 (2)<^ 1/2 (2)] 

-LX?/ 2 (1)0 P - 1/2 (1) + X A / 2 (l)^i/ 2 (l)][x? /2 (2)< /2 (2) + X A /2 (2)0 A / 2 (2)]} . (0.4) 

At this stage we use MATHEMATICA |||. We introduce Eqs. ([DD and the equivalent 



for the flavor parts in (p.4j) . We get a huge expression with 338 terms depending now on 
the quantum numbers of the quarks. In the matrix element of an operator O we then get 
338 2 = 114244 terms of the form 

(siS 2 S3S4S5S 6 rir 2 r3r4r5r 6 |0|s / 1 4s3 s 4 s 5 s 6 r i r 2 r 3 r 4' r 5' r 6) > (0-5) 
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where Si and 7* (i = 1, . . . , 6) stand for the spin and isospin projection of the i th quark. Note 
that the normal order of particles is implied. Now let us choose O = B\ ■ c?3 A{ • A3 P-$l -, 
which contains the permutation P36. Then we have 

(SiS 2 S Z S A S^S & T l T 2 T ?> T i T b T & \a l ■ a 3 \{ ■ A$P^ l*!^^^^ 7 ! 7 ^) 
= (S1S2S3S4S5S6T! T 2 T 3 T 4 T 5 T e | O x ■ O z \[ ■ ^S^s'^s'^T^TqT'^T^) 

= (*iW 3 |c?i ■ a 3 X{ ■ A 3 V 1 4^6> <J^<<»<<<<" 

= <«i*s|*i ■ WAHtiTiM ■ %\«) S|<<<<<<<' ■ (°- 6 ) 

This shows how a six-body matrix element can be reduced to the calculation of two-body 
matrix elements. The necessary nonzero two-body matrix elements are 

(TTI*i-*2|TT) = <UI*i-*2|U> = i, 

(Tl |*i • * 2 | Tl) = (IT l^i-^l IT) = -1, 
(TII*i-* 2 |IT) = (ITI*i-*2|TI) = 2, 

(uu\\{ ■ \ f 2 \uu) = (dd\\{ ■ \{\dd) = 4/3 , 

(ud\\{ ■ \ f 2 \ud) = (du\\{ ■ X f 2 \du) = -2/3 , 

(ud\X{ ■ \{\du) = (du\X{ ■ \{\ud) = 2 . (0.7) 

MATHEMATICA is then used to compute systematically the sum of the 114244 terms 



stemming from Eq. (|0.4j) . 

In Tables I & II all required six-body matrix elements obtained by this technique are 
listed. 
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TABLES 



TABLE I. Matrix elements (a\0\(3) of different operators O for (S,I) = (1,0). 



Q: 


NN 


NN 


AA 


NN 


AA 


CC 


P 


NN 


AA 


AA 


CC 


CC 


CC 


1 


972 





972 








972 


^36 


-12 


48 


12 


-144 


288 


-756 


A1.A2 


-2592 





-2592 








-648 


^3- A 6 

















-1296 


A l- A 2 -^36 


32 


-128 


-32 


384 


-768 


72 


\c \c p/o-c 
A 3- A 6 -^36 


-64 


256 


64 


96 


-192 


1152 


\c \c p/o-c 
A l- A 3 *36 


32 


-128 


-32 


384 


-768 


720 


A l- A 6 -^36 


32 


-128 


-32 


-48 


96 


720 


\c \c p/o-c 
A l- A 4 r m 


-16 


64 


16 


24 


-48 


1260 


0\.02 T\.Tl 


4860 





972 








108 


0"3-CT6 73 -T~6 


-900 


576 


1980 








1116 


o\.a 2 T1.T2 -P 36 


-444 


48 


12 


-720 


288 


588 


03 -ere T 3 .r 6 Pj[ 6 


708 


48 


1596 


240 


672 


-1092 


CTi.CT 3 ri.r 3 P 3 J 6 


132 


336 


12 


-720 


288 


-420 


cri.cr 6 ti.t 6 P 3 J 6 


132 


48 


12 


336 


-96 


-420 


ai. cr 4 ri.r 4 P 3 J 6 


36 


-144 


-36 


228 


288 


-1260 


01. 02 A{.A£ 


4536 





1296 








-18 


0-3.0-6 a^.a{ 


-864 


576 


1584 








1020 


0-1.0-2 A{.a£ p(q C 


-376 


64 


16 


-672 


384 


706 


0-3.0-6 x{.x{ p(q c 


784 


32 


1520 


216 


528 


-1024 


0-1.0-3 a{.a{ p(^ c 


104 


304 


16 


-672 


384 


-332 


o-i.o-g A{.A{ p/^ 


104 


64 


16 


340 


-200 


-332 


\ / \ / d/ ctc 
0-1.0-4 A1.A4 -T 36 


44 


-152 


-32 


278 


164 


-1197 


\/,0 \/,0 

o"i.o- 2 X J { .A^' 


-648 





648 








-252 


\/,0 \/,0 

0-3.0-6 A^' .A 6 ' 


72 





-792 








-192 


\/,0 \/,0 Tjfvc 

0-1-0-2 X{' ..X J 2 ' P 3 J 6 


136 


32 


8 


96 


192 


236 


\/,0 \/,0 D /crc 

0-3.0-6 A 3 ' .A 6 ' P 3 J 6 


152 


-32 


-152 


-48 


-288 


136 


\/,0 \/,0 D /crC 

0-1.0-3 A{' .A 3 ' Pfc 


-56 


-64 


8 


96 


192 


176 


0-1.0-6 AJ' .A 6 ' P& 


-56 


32 


8 


8 


-208 


176 


\/,0 \/,0 D /crc 

0-1.0-4 Af .A4' P 3 J 6 


16 


-16 


8 


-20 


-248 


126 


factor 


1 


V5 


1 


V5 


1 


1 


972 


972 


972 


972 


972 


972 
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TABLE II. 


Matrix elements 


(a\O\0) 


of different operators O for (S 


1,1) = (0,1). 




a 


NN 


NN 


AA 


NN 


AA 


CC 


a 


NN 


AA 


AA 


CC 


CC 


CC 


i 


972 





972 








972 


pf ac 

-^36 


-12 


48 


12 


-144 


288 


-756 


A1.A2 


-2592 





-2592 








-648 


^3-^6 

















-1296 


\c \c pfcc 
A 1' A 2 -^36 


32 


-128 


-32 


384 


-768 


72 


\c \c p/o-c 
A 3- A 6 ^36 


-64 


256 


64 


96 


-192 


1152 


A l- A 3 -^36 


32 


-128 


-32 


384 


-768 


720 


A l- A 6 ^36 


32 


-128 


-32 


-48 


96 


720 


A l- A 4 -^36 


-16 


64 


16 


24 


-48 


1260 


0-1.(72 T1.T2 


4860 





972 








108 


o"3-cr 6 r 3 .r 6 


-900 


576 


1980 
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-444 
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12 


-720 


288 


588 
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708 
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672 


-1092 
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-420 
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48 
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-376 


64 


16 
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814 


a 3 .(TQ X f 3 .X f 6 P(£ c 
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32 
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232 
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(71.(73 A{.A{ Pg 7 ^ 
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16 


-672 
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-260 


(71.(76 A{.A 7 Pg 7 ^ 
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16 
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-248 


-260 
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-168 
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03.(76 A3 .A 6 
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96 
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20 


-328 


210 


factor 


1 

972 


V5 
972 


1 
972 


V5 
972 


1 
972 


1 

972 



23 






0.2 



-0.4 



-0.6 



-0.8 



-1.2 



0.5 



FIGURES 



k(fm" 



NN 

NN+AA 
NN+AA+CC 



FIG. 1. Si NN scattering phase shift as a function of k. The solid line shows the result for 
the NN channel only, the dotted line for the iViV+AA and the dashed line for the iViV+AA+CC 
coupled channels. 



24 



0.2 



NN 

NN+AA 
NN+AA+CC 



"D 



-0.4 



-0.6 



-0.8 



-1.2 



0.5 



k(fm" 



FIG. 2. Same as Fig. 1 but for the Sq partial wave. 
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FIG. 3. The relative wave function of Eq. ( |5.3| ) for the S± partial wave for k 
obtained in one channel (solid line) and three channels (dashed line) calculations. 
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FIG. 4. So an d Si NN scattering phase shifts as a function of the laboratory energy E[ a ^. 
The solid and dotted lines show the result corresponding to the GBE model and the dashed and 
dot-dashed lines that of the OGE model (see Ref. pl|). 
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